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Abstract 

We consider an N = 1 SU(N C ) SUSY gauge theory with Nf > N c matter 
multiplets transforming in the fundamental and antifundamental representations of 
the gauge group. Using the Konishi anomaly and a non- anomalous conservation 
law, we derive a system of partial differential equations that determine the low 
energy effective superpotential as a function of the mesonic and baryonic vacuum 
expectation values. We apply the formalism to the cases of Nf = N c and Nf = N c +1 
where the equations are easily integrated and recover the known results. We further 
apply the formalism to derive a system of partial differential equations to determine 
the low energy effective superpotential for the Seiberg dual theories. Finally we 
briefly discuss the associated matrix models via the Dijkgraaf-Vafa conjecture. 



1 Introduction 



Quantum field theories with Af = 1 supersymmetry were studied in great detail in the 
mid 90's, for reviews see eg. [TJ Ej. One of the problems addressed was that of the 
dynamically generated superpotential. Seiberg j3j showed that one can go a long way in 
determining the exact form of the superpotential for various theories using symmetries 
and holomorphy. 

Recently a new technique for evaluating low energy effective superpotentials has been 
developed by Dijkgraaf and Vafa [UEIIS1- They conjectured that the exact low energy 
effective superpotential for an Af = 1 SUSY gauge theory with an adjoint matter field 
could be constructed from an associated matrix model. The action of the matrix model 
consists of the superpotential appearing in the definition of the gauge theory and the 
variables of the matrix model are the superfields themselves. The conjecture then relates 
the planar diagrammatric contribution to the free energy of the matrix model to the exact 
low energy effective superpotential of the corresponding gauge theory. 

Three different proofs of the conjecture have been given, |Hl El- [Zj treated the 
theory as an Af = 2 theory broken to Af = 1 by the presence of a superpotential and 
then analyzed the gauge theory using the Seiberg- Witten formalism [TU1 ITT] . The low 
energy effective superpotential was constructed and shown to be the same as arising 
from the associated matrix model via the Dijkgraaf- Vafa conjecture. In jS] the gauge 
theory was analyzed perturbatively with the result that the diagrams relevant to the low 
energy effective superpotential were indeed planar and corresponded directly to planar 
graphs of the associated matrix model. In |l] a generalized Konishi anomaly equation was 
constructed which takes the exact form of a matrix model loop equation. The solution of 
the anomaly equations was shown to correspond to the solution of the associated matrix 
model loop equations in exactly the way conjectured by Dijkgraaf and Vafa. 

The purpose of this paper is to analyze the Af = 1 SUSY gauge theory with funda- 
mental matter only along the lines of the Konishi anomaly analysis of [9 . This particular 
topic has been studied in [T3J CH E3 EH HZj. Specifically we shall be interested in 
theories with the number of flavors satisfying Nf > N c and vanishing superpotential. For 
such theories the low energy degrees of freedom consists of mesons and baryons. For 
Nf = N c and Nf = N c + 1 the exact results are known, JS]- I n the former case there is no 
dynamically generated superpotential but the classical constraint equation on the meson 
and baryon vacuum expectation values is modified. In the latter case there is a dynam- 
ically generated superpotential which serves to enforce the classical constraint equations 
on the meson and baryon vevs. 

The remaining cases Nf > N c + 2 are less well understood. For 3N C < Nf the theory 
is infrared free while for 3/2N c < Nf < 3N C the theory is believed to be a conformal 
field theory. In the range N c + 1 < Nf < 3/2N c the theory can no longer be conformal 
as such an assumption would violate unitarity. Instead Seiberg has conjectured [TH] that 
the theory has a dual description in terms of an SU(Nf — N c ) gauge theory with Nf 
fundamental matter multiplets along with a set of scalar multiplets. This conjecture has 
passed various consistency checks including showing that both theories have the same 
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global symmetries, have matching quantum numbers for the gauge invariant operators, 
and have matching anomalies. 

By deriving the Konishi anomaly equation, along with a non-anomalous conservation 
law, we construct a set of partial differential equations for the low energy effective su- 
perpotential W e ff as a function of the meson and baryon operators for Nf > iV c . The 
elementary superpotential of the theory consists of terms giving vevs to the mesons and 
baryons of the theory enforced by introducing Lagrange multiplier fields. The Konishi 
anomaly and conservation equations are enough to solve for the Lagrange mulipliers in 
terms of the meson and baryon operators. These solutions are then used to construct a 
set of partial differential equations for W e ff determining it up to a function of the glueball 
superfield. 

In sections 2 and 3 we present the derivations and results for the Nf = N c and 
Nf = N c + 1 cases respectively. In these cases we can solve the equations completely 
and are able to recover the results of 18J. In sections 4 and 5 we derive the Konishi 
anomaly and conservation equations for the theory and its Seiberg dual respectively for 
Nf > N c + 2. In this case the equations become considerably more complicated and we 
are unable to solve them. Nevertheless the set of equations could in principle be used to 
provide a proof of Seiberg's conjectured duality. In section 5 we consider the associated 
matrix model and derive its loop equations. 

2 N f = N c 

The situation with Nf = N c = N is somewhat simpler, so we begin with it. As our 
superpotential we shall simply take 

W = m f f (M f f-Mt)+b(B-B)+b(B-B) (1) 
where we have defined the meson and baryon operators 
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The fields m ^ 6, and b are Lagrange multipliers that enforce the conditions that the 
vacuum expectation values of the meson and baryon fields are given respectively by M*~, 

B and B. In the context of matrix models and the Dijkgraaf-Vafa conjecture, this kind of 
superpotential was introduced in t 20j for the mesons and generalized to include baryons 
in |2U. This theory has an SU{N f ) Q x SU{N f ) Q x U(l) R x n£i U(l) i>Q x n£ U{1\ Q 
symmetry 1 with the transformation properties of the fields given in the table below. 

1 Note that the U(l) flavor symmetries are not all independent because of the SU(Nf) symmetries. 
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If there is no superpotential then the flavor symmetries would classically give rise to the 
conservation equations D 2 (Q^e v 'Q?) = and D 2 (Q^e v Q ?,) = 0. The U(l) symmetries 
are however anomalous and are modified at one- loop |221 E9 to 

D 2 (Q\,e v Q f ) = S5 f fl (5) 
respectively, where the glueball superfield S is defined as 



D 2 (Q f e v Q n ) = S5 f ~ (6) 



S := -±-TrW a W a . (7) 

OZ7T z 

Reinstating the superpotential (JIJ, the anomaly equations pick up some extra terms from 
the equations of motion for Qf, and Qf'^ and become 

D 2 (Q),e v Q } ) = S5 f fl + M f f m f f , + bB5 f fl 

D 2 (Q f e v Q^) = SSj' +m f f 'M f f + bB8 f ~ (8) 

respectively. Taking expectation values in a supersymmetric vacuum will remove the 
left-hand-sides of these equations (or equivalently working in the chiral ring of operators 
modded out by D exact operators). We finally arrive at 

(M^)m f 'f, = ((S) - b(B))5 f f , 

m r f {Mt) = ((S)-~b(B))5j'. (9) 

We can obtain another set of equations similar to the above set by noting that there 
is another, gauge invariant, conserved current. Define the pair of operators 

R ai — 1 ai-a N Qh . . . nfN — p>(n- l \ ai 

n fi (jV-1)! 6 /i-/ivVa 2 Vajy-^lV >h 

A fl 1 A 

r •— e . . .n a ~ N — RfcHV 1 r ioI 

D at ■ (N-l)l. l N ^ fx 1 )ai ' 1 J 

Under SU(N C ) these transform exactly as Q°~ and Q{, respectively and therefore in the 

classical conservation laws above we may replace Q°~ by Bf and Q } a by B a to obtain new 
conservation laws 2 . Unlike the cases above however, these new conservation laws are not 



2 Such generalizations have also been considered in ^3] 
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anomalous (as there is no propagator for Q< and Q or for and Q) . Including the terms 
arising from the superpotential we obtain 

D 2 (Q\,e v B ) = m f f B + bBB(M~ 1 ) f f 

D 2 (B f e v Q f ~' r ) = m f ' f B + bBB(M- 1 ) f ' f (11) 

which may equivalently be written as 

* f 

D 2 {Q\,e v B ) = m f f B + bdetM(M- l ) f f 

D 2 (B f e v Q f '"<) = m ! f B + bdetM(M- l ) ! f (12) 

where now the determinant is taken over the flavor indices. Taking expectation values as 
before removes the left-hand-side. 

It was shown in [5] that an expectation value of a product of gauge invariant, chiral 
superfields is equal to the product of the individual expectation values. Applying this to 
the expectation values of ()12j) results in the equations 

m f f (B) +bdet(M)((M)- 1 ) I f = 

m f f (B) +bdet(M)((M)- 1 ) f f = 0. (13) 

We must be careful in obtaining this result. We could have taken the expectation values 
of these same equations in the form of (fTTj) . Applying the factorization property of [9 J 
would have resulted (after some trivial manipulations) in the same equations as in (|9]) 
althought without the S dependent term. Taken together these two sets of equations 
would have implied that S must vanish, which is not correct. 

Our goal now is to use these equations in order to solve for the low energy effective 
action W e ff(S, M, B, B). To do this we note that 3 
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Therefore if we can solve the above equations for the Lagrange multiplier fields in terms of 
the variables appearing in the effective action, then we can integrate to find the effective 
action up to an S dependent term. 

To begin we note that (0) implies that 

bB = bB (15) 



3 To simplify the notation we shall not write the expectation value signs (• ■ •) from now on. 
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and therefore that 

B = ~bi, B = b£ (16) 
where we have introduced a new variable £. We may therefore rewrite as 

M f f =(S-!j.){m- 1 ) f r (17) 

We have at this stage reduced the problem to solving for £ in terms of the meson and 
baryon vevs. 

To find an equation for £ we need only eliminate the Lagrange multipliers m f and b 
(or b) from (|13|) using the expressions in (J 17)) and (|16|) . This results in a linear equation 
for £ with solution 

det M — BB 

« = ~ s • (18) 

Substituting for £ in (|Td^ and (fTT)) we find 

6 = -S- 6 



det M -BB 

6 = -S 

det M - BB 

m* = S d6tM A M-% (19) 
/ detM-fifi V Jf K J 

Finally we can insert these expressions back into the derivative equations (|T3jl for W e ff 
and integrate. This results in 

W«, = -S* f* M K - BB )+m (20) 

where f{S) remains to be fixed and we have put in the scale A in order to make the 
argument of the logarithm dimensionless. f(S) is just the Veneziano-Yankielowicz con- 
tribution 

W V y(S) = (N c - N f )S(- In ^ + l) (21) 

which vanishes in this case as iV c = iV^. Therefore setting f(S) to zero and minimizing 
W e ff with respect to S results in the constraint 

det M — BB = A 2N (22) 

and moreover a vanishing superpotential. This is in agreement with the results of [THj . 
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3 N f = N c + 1 

The methodology in the case Nf = N c + 1 is identical to the previous section, although the 
details are slightly more involved. As before the superpotential that we consider consists 
of Lagrange multiplier terms that fix the vevs of the gauge invariant variables, i.e., 

W = m f f (M f f - M f f ) + b f (B f - B f ) + b f ~(B - B f ~) (23) 
where the meson operator was defined before, the baryon operators are defined as 

Bf ■■= jtflh--l„<? l ~°"Qk---Ql"„ (24) 

and the fields rrij, b^, and bj are Lagrange mulitpliers. 

Proceeding exactly as before, we must derive the "loop" equations analogous to those 
found in Q and (fI3|). In the case of © the steps are virtually identical, so we shall forego 
the details and simply record the result. After taking expectation values we find 

{M f -)m ! f , = ({S)-V>{B g ))6 f fl + bf{B f ) 

m^Ml) = ((S)-b~ g (B))8l' + b f (B ). (26) 

The second set of equations follows analogously to those found in (|T3*j) . Namely we 
note that Ql transforms under SU(N C ) in the same way as the operator 

B? ■= ^f'^e aa2 ... aN Ql ■ ■ ■ Q% (27) 

t v~' S ' f 

and therefore for vanishing superpotential we may define a new conserved current (Q\ f e B ), 
which is not anomalous. In a completely analogous way we note that QJ transforms under 
color in the same way as 

Bff> ■= ^ff'h-fn^-^Q^ ■ ■ ■ Qt (28) 

and therefore for vanishing superpotential one has the conserved current (Bff>e v Q^^). 
Reinstating the superpotential and taking expectation values results in the equations 

(N - iy.(m f f (B ) - m? f (B)) + ^e / ,.. /jv _ l/s e^-^-^'(M^) • • • (M^> = 
(N - l)\(mf f (B f/ ) - m } f ,{B f )) + 6 i e^-^-^e /l ... /w _ l// ,(M|) • • • (M^> = (29) 

where as before we have used the factorization property of jH]- 
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Now we can proceed to solve these equations. Recall that the goal is to construct the 
low energy effective action W e ff(S, M, Bf, w) where W e ff satisfies the partial differential 
equations 



dW, 



dW eff 
dB f 

dBf 



— — —rrJ 



f 

-bf. (30) 



We can use the equations (j2Tlj) and (}2T?j) to solve for rrif, bf and bj in terms of S, M p Bf 

and Bf , and then integrate to find W e ff up to an S dependent term. 
To proceed let's first rewrite (|26J) in the form 

Aft = {S -b 9 Bg){m- 1 )j + b f B r {m- 1 ) f f ' 

M f f = (S -bgB^m^Yj + bfim-^B (31) 

where we have dropped the expectation value signs and will for the remainder of this 
section. Taking the trace of this equation results in the identity 

~ a,/ 

b f B f = b?B . (32) 



Inserting this back into (j31|) leads to 



b f Bf^m-y'- =b f ~B (m _ %. (33) 



This equation is enough to imply that 



Bf = bpnk 

B = mf f b f £ (34) 

where we have introduced the unknown £. One further result that we need from (J31)) can 

s./ 

be obtained by contracting it with either BfOiB . After some simplification and use of 
the relations we find 

(m -1 )0 = ~M } f B 

Bf(m- l ) f f ~ = ^B f M f f . (35) 
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Using (J54*j) and (J33j) we may rewrite (|3~T]) as 



= (5 - ^^)(M-y f - ifiB, (36) 

Combined with ()34|) we have managed to express the Lagrange multipliers in terms of the 
desired set of variables up to the unknown £. 

To solve for £ we must now bring in the equations (|29)l. Eliminating mj and 6 9 or 

bg via the relations ()36|) and respectively results after a bit of algebra in a linear 
equation for £ with solution 

det M - B q M%B , s 

£ = g 2 9 9 ■ (37) 

This finally allows us to express the Lagrange multipliers in terms of the desired variables. 
We can now integrate the equations (replacing the operators M etc. by their vevs) for 
W eff (dlj) to find 

/det M-B q M% B\ 
Weff = -8hx( A J 9 ) + f(S) (38) 

up to the S dependent function f(S). f(S) is given by the Veneziano-Yankielowicz term 
(|2T|) for Nf = N c + 1. Substituting and integrating out 5* we find for the low energy 
effective superpotential 

Weff = -^=5 (det M - B g M%&) (39) 
in agreement with the results of [TH] . 



4 AT/ > iV c + 1 

The results obtained thus far are not new, only the methods to obtain them are somewhat 
novel. The more interesting cases in which to apply these methods are when the number 
of flavors is even larger, i.e., when Nf > N c + 1. In this case one can hope to address 
questions about confinement, chiral symmetry breaking, and Seiberg duality [T^. While 
the methods used above are straightforward to apply to the cases of Nf > N c + 1, the 
equations that one obtains are more complicated and difficult to solve. We shall outline 
the construction in this section and leave its possible solution to future work. 

As before we consider a superpotential which consists solely of Lagrange multiplier 
terms enforcing the meson and baryon vevs, 

W = m f f (M - M f f) + // • / m7V a " B h ... f -) + h J: ... f il/ ' U - B h ~**) (40) 
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where we have defined N := Nf — N c = Nf — N and the baryonic operators 

Bt t ■= —e f f e ai "' aN O gi ■ ■ ■ 9N (41) 



AT! 

R ■= e h~f$9i-9N e 0~ 1 ---0- N (42) 



The remaining fields in the superpotential, rrif, b^ 1 ""^ and are the Lagrange mul- 

tipliers. 

As before we first construct the Konishi anomaly equations which are given by 

(M f f )m f fl = ((S) - b*"f*(B fl ... ff .))8 f fl + M^*-if'(B fl ... f ,_ if ) 

mfrM'j) = ((5)-6 /i ... 4 (5 r '^))4' + 6 /i ... / ^ i/ ,^ /l "^- l/ ). (43) 

The other pair of equations needed to solve for the low energy effective action follow in 
direct analogy to before. We define the operators 

r •— e h-fa +1 92-9N e O a - 2 ---O aN 

R a — — f , / p <H*2-aN rite . . . fpN (AA\ 

We note that these transform under the SU(N C ) symmetry exactly as Ql and QJ respec- 
tively and therefore for vanishing superpotential one has the exact conservation laws 



D\Q)e v B '"^) = 
D\B fl ... fA+i e v Q") = 0. (45) 

Putting the superpotential (}4T)j) back in and taking expectation values results in the equa- 
tions 

(N l)\m f f ((B )& r -Y,{B ><f) 

i=i 

AT 

(N l)\m f f ((B gi ... 9 ^ g -j:(B 

K) 

i=i 



3-/i"'/]v+i 



where we have used the factorization property of gauge invariant, chiral operators in 
writing these equations. 

The procedure for solving these equations for the Lagrange multipliers in terms of 
the meson and baryon expectation values is in principle straightforward. We begin by 
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multiplying (|46|) by (N+l) M operators and contracting with the free upper/lower indices 
respectively. This produces a det M operator in the second term. Next we use (f4*3j) to 

eliminate the mj- Lagrange multiplier in (jUj) . After some algebra we arrive at 

AT N \ 



b 9v " 9 * = -^-^ — H (S+ -2)(b- B))(MB) 9i - 9 n 

N & * - - \ 

- (-£ - l) • B)? {M^YfMV (MB)Bi-~Bi-ifBi+i~Bii j (47) 

j=l ) 

where we have introduced some notation to simplify the expression, 

(ft.fl) := ^-^4,.. 9J , (48) 

with the other quantities defined similarly. 

At this stage we can eliminate one Lagrange multiplier, b 91 " 9 ^ say, to find a linear 
equation for bg v ..g- in terms of the meson and baryon expectation values. Unlike the 
previous cases however the resulting equation is considerably more difficult to solve. The 
key difference is that the last terms on the right-hand-sides of the equations in (j4*7|) were 
of the same form as the first terms on the right-hand-sides. In that case solving the 
equations was trivial because they could be reduced in the end to solving for a single 
variable. In this case no such reduction (as far as we know) is possible making solving 
the equations more involved. 

One interesting point to note however is that the solution for the Lagrange multipliers 
is necessarily linear in the gluino condensate S. This is easy to see by simply scaling all 
Lagrange multipliers in (J47|) by S and noting that S then factors out of the equation. The 

same holds true in the equations for in (|43|) . As a consequence, after adding in the 
Veneziano-Yankielowicz contribution, the low energy effective superpotential will take the 
form 

W eff = (N c - N f )S(- In A + i) + S f(M, B, B) (49) 
where f(M, B, B) denotes the contribution coming from solving the equations above. 



5 Seiberg dual theories 

When the number of flavors satisfies N c + 2 < Nf < (3/2) iV c , the moduli space of 
the gauge theory discussed in the previous section is known to be singular at the point 
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M-j = Bf v ..f. = B^'fa = 0. It was conjectured by Seiberg that the singularity is due 
to the presence of particles becoming massless at this point in moduli space. Seiberg 
conjectured a dual IR description of this theory which moreover is valid at the origin of 
moduli space. This theory consists of Nf chiral matter fields q$ and q* , a set of chiral 
scalar fields M~, and superpotential W m = M~ qfqf where the gauge group is SU(Nf—N c ) 
and the gauge indices have been suppressed. 

There are several non-trivial consistency checks that are discussed in detail in [TS] that 
this dual description must pass in order to have a chance to describe the same IR physics 
as the theory discussed in the previous section. First note that the global symmetries are 
the same. Secondly we note that the massless degrees of freedom away from the point 
= are the same. To see this let us first define the dual mesonic and baryonic 
operators 

N f f := qfq f 

B (m) := e<*- a *qf{---qf* (50) 

where the (m) super or subscript denotes that this is the dual baryonic operator. Because 
of the superpotential the dual meson is massive away from the point M~ = 0. The baryonic 
operators of the dual theories have identical transformation properties, identifying them 
then the dual theories in the end are described by the same set of low energy degrees 
of freedom. Finally one can check that the 't Hooft anomaly matching conditions are 
satisfied. 

To go even farther one would like to compute the low energy effective superpotential on 
both sides of the duality and show that they are the same. While we will be unable to do 
this we nevertheless think that it will be useful to use the analysis described above to find 
a set of equations whose solution would provide the low energy effective superpotential. 
The analysis is essentially identical to above, the only new features here are that in this 
case there is already a superpotential present, aside from the Lagrange multiplier terms 
that we will add, and also that this time around we have an extra set of scalar chiral 
matter fields to deal with which leads to more equations. 

The superpotential that we will take is 

W m = M f f N? + rrJ f {M f f - M f f ) + l/ 1 "** {B^}. f , - B { ™}. ff . ) 

+ h^w'* - Bfc'*) (51) 

where as before the Lagrange multipliers rrif, b^ 1 '"^ and enforce the vev con- 

straints. This form is valid provided that the vev of the scalar fields is non-zero, 
which is all that we shall concern ourselves with here. The analysis now proceeds as 
before. In complete analogy to the previous sections (essentially one just replaces Q's 
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with g's) we find the set of Konishi anomaly equations (dropping expectation value signs 
as before) 



M f f N f fl = sdj, ■ .Y//-A :/ //;.; , :; A . r 



and the set of conservation equations 



MB^f. f = -N\b f f fM-'-N^- 1 . (53) 



Finally we can get another set of equations by recalling that M - is a chiral, gauge invariant 



scalar field. Taking the expectation value of its equation of motion leads to 

J f + Nj = 0. (54) 



The strategy before was to solve for the Lagrange multipliers as functions of the low 

energy degrees of freedom M-~, Bj x ,,,^ and B^ . This leads to a set of equations in 
complete analogy to (j3~Uj) which can in principle be integrated to find W e ff up to an S- 
dependent term. As an example, we can easily use the equations above to construct an 
equation for fr^.../.- Multipy the second equation of (|53j) by an M on every free index. 

Then use flSlj) and (J52J) to replace the N dependence in terms of the desired variables. 
We find in the end 



= M f }, ■ ■ ■ M f f B^ , 
fi fx n * N 



N 

f N— 1 707/7 TO \ JS- 



+ N\b, } . ?(S8 f l-N{b-B im) )p)--.(S5 T *-'-N(b-B (55) 

N 1 N J 1 J 1 N— 1 N—l 

i.e., a nonlinear equation for b in contrast to the linear equation that we obtained previ- 
ously (pT7j) for the IR dual theory. 

It is interesting to note that some information on the S-dependence of W e ff can be 
obtained without solving for W e ff explicitly. If we scale the Lagrange multipliers b and 
b as well as the mesons by S, then all S dependence divides out of (|53|) as well as 
from Consequently we know that W e ff must take the form, up to an S'-dependent 

constant of integration, W e ff = Sg(M/ S, B, B) for some function g. Putting in the 
Veneziano-Yankielowicz contribution we find 

W e ff = -N C S( - In + l) + Sg(M/S, B, B). (56) 
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6 Matrix model approach 



In light of the conjecture of Dijkgraaf and Vafa, it is natural to ask if the results ob- 
tained earlier can be obtained from a matrix model. This question has been discussed 
by various authors [2U EJ 123 EH] • Recall that the Dijkgraaf- Vafa conjecture relates the 
planar diagrammatric contribution to the matrix model free energy to the low energy 
effective superpotential of the gauge theory. Moreover the matrix model action is given 
by the gauge theory superpotential expressed in terms of the elementary superfields. The 
question becomes somewhat more subtle when baryons are involved, for two reasons. The 
first is that normally one extracts planar diagrams from the matrix model by taking a 
large N limit, that is, the planar diagrammatic contribution dominates in this limit with 
all planar diagrams scaling with N in the same way. When baryons are included, the 
diagrams that one would like to call "planar" 4 no longer all scale with N in the same way, 
making their extraction complicated. A second reason that baryons are more troublesome 
here is that the previous discussion assumed that the number of flavors Nf was also taken 
large, i.e. Nf = N + k where k is fixed, but it has been shown that when fundamental 
matter is included, only planar diagrams with at most a single boundary are relevant 
to the low energy effective superpotential. Taking Nf large will not help to isolate the 
planar diagrams with different numbers of boundaries. Nevertheless, it has been shown in 
|2T] that the generalized planar diagrams (in the gauge theory) are the relevant ones that 
give rise in the end to the low energy effective superpotential. If it is also true that the 
generalized planar diagrams from the associated matrix model are relevant to computing 
W e ff, then extracting them remains another question. 

In this section we derive a set of matrix model loop equations for arbitrary iV c and 
Nf > N c . These equations have identical form to those constructed in previous sections 
provided that we identify the glueball field S with g m /N where g m is the matrix model 
coupling constant, except that we don't have the factorization property that followed in 
the gauge theory case from chirality. Normally on the matrix model side this property 
follows from the large N limit, namely in the large N limit the disconnected diagrams 
dominate the correlator giving rise to factorization. If one can make sense of the large N 
limit in this case, then this might be good enough to derive the factored form of the loop 
equations thereby giving identical results to the gauge theory analysis above. However 
the naive large N limit does not produce sensible results as the different terms in the loop 
equation scale with N differently. 

Deriving the matrix model loop equations is a straightforward matter. As an example 
we shall consider the N c = Nf case. Our matrix model is given by the integral 

e~ N2 F = J [dQ] {dQ]e-W^ s (57) 

with action S given by the superpotential (JTJ) considered previously in the Nf = N c case. 
Here the matrix model variables Ql and QJ are N x N where N c = Nf = N. To derive 

4 Some generalization of the definition of planar is needed here due to the presence of the e tensor in 
the definition of a baryon, see eg. [21] • 
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the loop equations we consider first the identity 

= J \dQ][dQ\-^ l (Qf; C -W*»>*). (58) 

Expanding this out results in the equation 

{Mt)m f f , = ((S) - b(B))5 f f , (59) 

where we have identified S = g m /N. Replacing the explicit Q variables by Q variables 

in ()58|) results in another loop equation of the same form in terms of the b and B fields. 
These loop equations are in fact identical to those in the gauge theory case. 
The second set of loop equations is obtained by considering the identity 

= j [dQ] [dQ]^f {^f^J h -- jN Ql ■■■Qf N e" (A ^ )s ) . (60) 
Expanding out the derivative results in 

m f f (B) + 6(det M^M" 1 )^) = 0. (61) 

Exchanging the Q and Q variables results in a similar type loop equation involving b and 

B. The important point to note here is that the second term is an expectation value of a 
product of operators. To reproduce the gauge theory results we must replace this with a 
product of expectation values of individual mesons. 

As we said before, this factorization is usually accomplished by taking a large N limit, 
however the different terms in the above loop equations scale differently with N making 
this limit less than helpful. Perhaps there is a double scaling limit that one could take 
that preserves the form of the loop equations while giving rise to factorization. 
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